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Abstract 

We have studied a Fermi system with attractive U(V)-symmetric interaction at the finite temperatures by the quantum 
field renormalization group (RG) method. The RG functions have been calculated in the framework of dimensional 
regularization and minimal subtraction scheme up to five loops. It has been found that for r > 4 the RG flux leaves 
the system’s stability region - the system undergoes a first order phase transition. To estimate the temperature of 
the transition to superconducting or superfluid phase the RG analysis for composite operators has been performed 
using three-loops approximation. As the result this analysis shows that for 3 D systems estimated phase transition 
temperature is higher then well known theoretical estimations based on continuous phase transition formalism. 
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Introduction. 


The investigation of quantum Fermi systems and the phase transitions in these systems are the problems of perma¬ 
nent interest. To describe the quantum equilibrium Fermi system we use the temperature Green functions formalism, 
quantum field theory methods and the renormalization group approach. The analysis is based on the microscopic 
model with local attractive interaction of the “density-density” type ID El El- The model’s field action has the form 


S * - ^ 


Wa(d, 


A 

2m 




(1) 


where i// a , if/„ describe the fermion fields at the finite temperature T, these fields are complex-conjugate elements of 
the Grassmann algebra and a — 1,..., r, r is the number of spin degrees of freedom, A is Laplace operator; m is a 
mass of the particles; yu is the system’s chemical potential; A = An\a s \/m is positive coupling constant and a s is the 
scattering amplitude for interparticle 3/J-scattering; t is the “imaginary” time and t e [0,/? = 1/T]. All the necessary 
integrations and summations in formula 0 and similar expressions below are implied. It is also necessary to impose 
the antiperiodic boundary conditions with respect to “imaginary” time on the fermion fields. 


<Mp.O) = -tMp,j6), ^(p,O) = -*t(p,0. (2) 

In the r = 2 case, this action 0 corresponds to the Bardeen-Cooper-Schrieffer theory and a =j\ i are the two 
possible spin projections. The theory describes low temperature superconductivity in electron systems. We will 
consider the case of arbitrary even values r. It can be corresponded to systems of high spin fermions investigated 
recently 00®EiEii5i[ini[ni or with the electrons in solids which have a sublattice index and/or index corresponding 
to a degenerations of zone structure m. 
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Usually in the system under consideration the phase transition temperature is determined by the appearance of 
an anomalous solution of the Dyson equation fT|, and the order parameter of a superconducting phase transition is 
given by means of the composite operators yp a t// y j and . To investigate this model using renormalization group 

method the action is transformed by introducing the new boson fields \,X 1 USES). The action of the form 


= ifrUd, + Sp)l// a + tr^ 1 ' - ^tlxoryll>\ - Tjlf'axiy'/' 


ayYy 


( 3 ) 


was considered, where e p = p 2 /(2 m) - /i. 

It can be easily proven that the integration exp(-.S', //A ) over the fields x,x' leads to exp(-S^). The new fields are 
complex skew-symmetric matrices because the fields ip, t// ; are Grassmann variables. The Schwinger equations 

(Xly + WWy) = 0 , 

C Yya 4\p Q <\Jjy) — 0 . 

show that the x,X 1 determine the order parameter of a phase transition. 

The integration of expi-.S' t , Ar ) over the fermion fields ip, <// leads to the new action for the boson fields x^X' 


s x = tr UU " trln 


~X f 



( 4 ) 


here u> s = jtT(2s + 1) are Matsubara frequencies, s e Z. Using the Taylor expansions for ln(l + ...) we can rewrite 
the action as 



here wave lines denote the field x\ the plain lines denote free (tppj propagators, cross corresponds to i// : , x* 
fields. 

To obtain the effective action in the infra red (IR) region we have to present ([5J in the form of a Ginzburg-Landau 
functional by expanding all diagrams in the external momenta p and frequencies. Then x\X fields can be considered 
as f-independent. As a result the effective action has the form 


s x = tr^fcop 2 + t 0 )x + ^ tT(xx f )ti(xx f ) + tfCurW)- 


g 02 . 


(6) 


The term with got coupling constant was included to obtain the multiplicatively renormalized theory. The parameters 
of the action c o and go 2 are positive and they can be calculated from the expressions 
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here D is dimension of space, oj s are Matsubara frequencies. The integration over k is performed in a narrow neigh¬ 
borhood of the Fermi surface |e k - /r| <6. The parameter 6 can be similar to Debye frequency a> D for the system of 
electrons in the solids or similar to Fermi energy s F for ultra cold atoms systems and 6 = (2/e j lr ’’s l -- ~ QA9e f for 
such systems Ifl3l . 

The IR behavior of the model was studied in mm. The renormalization group (RG) investigation in 
the framework of s = 4 - D expansion in one-loop approximation m and then in three-loop approximation 03 
establishes the absence of IR-stable fixed points for even values of r > 4. It was found that the stability criterion for 
action ([6]) (the condition for positive definiteness of an interaction) can be formulated as the inequality 


g 2 + rgi > 0, 
2 


( 9 ) 






for g 2 > 0. 

Moreover solutions of the RG equations for the invariant charges in one-loop approximation m show that the 
system loses the stability before the continuous phase transition occurs. It was supposed that a first-order phase tran¬ 
sition takes place here,/ and this phase transition can be considered as one of the possible reasons of high temperature 
superconductivity. 

Then the similar behavior was confirmed in lfl5l in the three-loop RG analysis of the 3D and 2D models. But it 
was found that the three-loop approximation is not sufficient to ensure an accurate calculation of the phase transition 
temperature. Therefore we have to develop our analysis up to five-loop calculations, that is the maximal order available 
now in the framework of s- expansion fl9l . 

In Sec.2 we describe the five-loop RG analysis of the model investigated with r > 4. According Ifl4l there is no 
IR stable fixed point in the framework of s expansion. Thus instead of seeking fixed points of the RG equation we 
restrict ourselves to analysis of phase trajectories. It was indicated in DU that the equations for the invariant charges 
can be constructed in the e-expansion form. 

To analyze the phase portrait of these equations in the physical space dimensions (e = 1, e = 2) we must resum 
the terms calculated, for instance, using the Borel resummation technique. Such a resummation requires knowing 
the higher orders asymptotics (HOA) of the e expansion. The HOA of the considered model was determined in jl5l 
using methods of instanton analysis da. The analysis and the results obtained are described briefly in Sec. 3. It is 
interesting to note that we have found several instantons with different matrix structures. These instantons are essential 
in a Borel resumming at different values of the charges gi, g 2 . 

In Sec.4 we resum and solve numerically the RG equations for invariant charges. It is confirmed that the invariant 
charges in the 3D model cross the boundary (|9]» of the stability domain of the action |[6}. As for the 2D model, it is 
found that five-loop approximation is not sufficient yet for the accurate description of the phase transition type. Our 
results show that the type depends on the initial value of the coupling constant gy). 

In Sec.5 the first order phase transition is studied in 3D and 2D model to find the real phase transition temperature. 
The additional terms (~ * 6 ) were introduced in action These terms are IR irrelevant for critical behavior, but are 
relevant for the first order phase transition description. They are renormalized as a composite operators in three-loop 
approximation. Their contributions to the state equation are Borel resummed and the phase transition temperatures 
estimated. 

1. Renormalization group analysis 

The renormalized action of the considered model is given by the expressions m 

Sr = z\ tr/(-A)* + Z T Z 2 rtrx f X + Z gl Z* (tr*/) 2 + Z g2 Z 4 trxxW- (10) 

This expression is obtained by the multiplicative renormalization 

goj -> gjM s Z gi , x ~^X z x , t 0 -> tZ t , ( 11 ) 

here the parameter M is a so-called renormalization mass; g\,g 2 are dimensionless renormalized coupling constants, 
index zero denotes bare parameters. In this paper we use the dimensional regularization, e-expansion and the minimal 
subtraction scheme (MS-scheme) ffiD. The bare parameters goj and tq are associated with the microscopic parameters 
(|8][7]i by the relations g 0j = 'g 0j /c t 0 = r 0 /c 0 and* -> */ g/cg. 

Let us introduce basic elements of the Feynman diagrammatic techniques for the model. In momentum represen¬ 
tation the free propagator has the form 


jy-h-fe 

rjlj2 _ <t‘2 


FL JU2 = 

k 2 


- T 0 


and \{6 hh 6 hh - 8 hh 8 ilh \ 


here <5,y is Kronecker symbol; k is the momentum or the wave vector (h = 1). The tensor W/ 1 / 2 is antisymmetric with 
respect to the transpositions of its indexes i\ <-> L and /j <-2 j 2 , and symmetric with respect to the transposition of 
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Figure 1: Tensor structures of the vertices and propagator: vertex I corresponds to gi, vertex II corresponds to g 2 , and vertex III corresponds to the 
propagator. 


the pairs (ij, i 2 ) <-> (j \, ji). One can write the tensor structures for the vertices gi and g 2 too, but we give only their 
graphical representation (FigJTJ , the vertices antisymmetrization is implied. 

In MS-scheme all renormalization constants have the form of the poles in s 


Ze 



e = (gj,T,x) 


here Z](g\ ,g 2 ) denotes the residue at the simple pole in e for the corresponding renormalization constant. Let us 

remark that the interaction (tixx^) must be included for the multiplicative renormalizability of the theory. It is easy 
to verify that the corresponding counterterms appear due to the renormalization of the theory starting with the simplest 
one-loop diagram. 

The RG-functions (the coefficients of the RG equation |fl6l ) are defined by the relations 

P gj = D M Z gi , y e = D m lnZ f , (12) 


here D M is the differential operator M8m at fixed bare parameters, p gj are beta-functions of charges g,, and the 
functions y e are anomalous dimensions for parameters e. In the MS-scheme the RG-functions are connected with the 
renormalization constants by the following expressions m 


Pgj = sM + Vs,)* ye = -gkdkZj,. 


( 13 ) 


Program “FORM” was used jT8l for the tensor structure calculations of the graphs. Tensor structure of the graph 
can be factorized, and the rest of the diagram is equivalent to diagrams of the scalar O 4 model, the values for these 
diagrams are taken from well known five loop calculations of the (9(«)-symmetric <I> 4 model |fT9l[20 'l. Finally, in the 
five-loop approximation (about 120000 diagrams), the RG-functions of the theory were calculated. The rescale of the 
charges g ,• —> g,/ \6n 2 was used. Results of our calculations were controlled for r = 2, r — 3. In these cases the model 
(10 1 is equivalent to the 0(2)- and 0( 6)-<J> 4 models with vector order parameter, respectively. 

The RG equation leads to the known equations for the invariant coupling constants 


£^=0 = gi> where f = ln ^2 • (14) 

The infra red (IR) regime f —> -oo is usually connected with fixed points (g* r g* 2 ) that are determined by the conditions 
Pgi(g*,g* 2 ) = 0 for all indices i. The fixed point is IR-stable, if the matrix = d g p g .(g*,g 2 ) is positively defined. 
However, it was found, in the one-loop approximation OH, that these points do not exist for r > 4. There is the 
IR-stable fixed point in the model at r = 2, this point describes the critical behavior of the superconducting phase 
transition in systems with 1 /2-spin fermions. 

We will not search for the possible fixed point in the five-loops approximations of the model considered, instead 
of this trajectories of the invariant charges will be studied in the next sections. 
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2. Instanton analysis 


Consider the equations ( fi~4| ) with /^-functions ( |48|49| >(see Appendix). After the scaling of the charges g, and the 
dynamical variable of the RG equation (; as g, —> eg , and £ —» %/s, we get eq. (fl4]> in the form 


deg, = -gi + ^ e N B) N \g 1 , g 2 ), i= 1,2 


N =0 


(15) 


gilf=0 _ 


Explicit expressions of the B^ ) (g\,g 2 ) can be obtained from (48 49 1 . In our case K = 4 (five-loop approximation). 
Equations (fT5j) can be solved in the form of e-expansion with the formally small parameter e. Similar to lfl4l we 


will consider numerical solution of the equations (15 1 . As usual, the e-expansion in the right hand side of equations 


(15 1 is asymptotic expansion with zero radius of convergence. Then the equations 0 must be resummed to obtain 
results at physical points e = 1 or e = 2. The resummation process requires knowledge about the asymptotic behavior 
B {N) (g\,g 2 ) at N —> oo. Such asymptotic behavior is called a higher-order asymptotic (HOA) and was investigated in 
lH5l in the model considered. 

Let us recall the main details of the analysis 03. The investigation of the asymptotic behavior of higher-order 
perturbation corrections proposed in IflTlI is based on the saddle-point expansion of the path integral (instanton ap¬ 
proach). Calculation method for the HOA for renormalization constants in MS scheme developed in Il22ll was used. 
Partially renormalized Green function was considered where subtractions of all the divergences in subgraphs up to 
order N - 1 are supposed. The coefficients of the expansion in the parameter s of the 2T-point Green function 

G 2k (s,xi,...,x 2k ) = W~ l J' DxDx ] x(x\)x(x 2 ) 1 .. ■ xi^k-Gx^kf e~ sR , (16) 

W= j DxDx^e~ SR . 

can be calculated in high orders (N —> oo) with the use of the saddle-point method in the integral representation ED 


(xi,...,x 2k ) 


(-1)^ f dsG 2k (s , xu..., x 2k ) 
2m J (-e)^ +1 
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(17) 


where y is a closed contour encircling the origin in the complex plane of e. As usual JT7J, we will find the HOA at 
t = 0 and D = 4. After the rescaling of the parameters gj —> g;/N, x —> V/Vy, x’ —> \Nx^ the variational equations 
for functional Sr + ln(-s) with respect to the field variables and s takes the form 


- At + ^-xtrxx f + ^yXX f X = o, 

- At 1 ' + Y~x f tr xx f + ^x f xx f = o, (18) 

f dx {- f 1 ( tr xx f )~ + ^ trxx f xx f j = - 1 - 

Similar to 03, the counterterms Z e - 1 in action (jT0|» are irrelevant for the calculation of the stationary points. For 
matrix fields x, A 1 we can assume without loss of generality the block-diagonal Phaff’s form containing of p = r/2 
blocks 


X = diag(sicr,..., SpCr), x 1 = - diag(s*cr,..., s*(r), cr = I” 


0 -1 

0 


(19) 


with some complex functions s/x). Any skew-symmetric matrix can be reduced to this form by some unitary trans¬ 
formations U(r). The equations (18 i and (19 1 yield the system of equations for sj(x) 


As,-(x) + sgi ^ \s k (x)\ 2 Si(x) + ^p|Si(x)| 2 s ; (x) = 0. 


( 20 ) 


k= 1 
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We seek ,v ; (x) in the form 


Si(x) = 


cti y 


ai G C, 


( 21 ) 


|x - x 0 | 2 + y 2 ’ 

similar to solutions of variational equation for the scalar O 4 model. Functions ,v,(x) depend on Xq and y - arbitrary 
parameters reflecting the translational and dilatation invariance of the theory. Then the Faddeev-Popov method was 
used similar to IfTTl . The method fixes the values of the free parameters for each realization of ,v,(x). Substituting pi) 
in (20 1 , we get the system of algebraic equations for constants a,- 


8 + egi Y ktlV- + -§-|(*;| 2 £h' = 0. 


( 22 ) 


k=l 


One can see that the stationary solutions may contain m - 0,..., p - 1 zero blocks with |a,| = 0 and n = /?,..., 1 
blocks with |a,| 2 = -1 6/(2neg\ + sg 2 ). and n + m = p. Phases of the complex numbers a, are not fixed, they are 
arbitrary parameters as well as xo and y. In addition to parameters y, xo and phase factors, there is also an invariance 
under unitary transformations U(r). Thus the number of zero modes is determined by r 2 - 2r + n + 5. This number 
influences the exponent of N in the HOA. 


Combining the instanton solutions (21 1 , the Pfaff’s form (19 1 and the third equation (181, we get stationary point 
in e parameter as 


sAn) = 


4n 


1 


(23) 


Similar to 
pole in s. 


3 2ng i + g 2 

, the beta-functions HOA can be obtained from the HOA for Green functions residue at the simple 


^\gugi) = constiN\N b - (-«)*(! + 0{n~ 1 )), 


(24) 


here consti - some constants not essential for future analysis, b n = (r 2 — 2 r + n + ll)/2 and a = max|a(n)|, a{n) = 

n 

-1 /E st (n). One can see from p3) , that a(n) depends on values of g;, therefore the largest of all a(n) gives the largest 
contribution to the HOA. Thus the perturbation series in the parameter s have zero radius of convergence in the theory 
with the action (fTO). For this reason, it is necessary to use some procedures of resummation e.g. the Borel method. 


3. Solution of the RG equations. 

3.1. Resummation of the RG equations. 

Let us recall the basic expressions for the Borel resummation t24l . We assume that there is a function Q(e) defined 
as a series on the parameter s 

Q(e) = Y bN Q {N) ’ (25) 

iV>0 


and the higher-order asymptotics of the series coefficients are determined by expression p4) . The Borel transform of 
the series i 


is given by the relations 


co 

Q(s) = f dt e-'t bo B(st), B(t) = Y B (N) t N , B (N) = 
n n> o 


Q 


. m 


r(N + b 0 + iy 


(26) 


where bo is an arbitrary parameter. The known asymptotic expansi on p4) together with several assumptions about 
the analytic properties of B(t ) allow one to resum series (25 i using (|26|) and to obtain a more precise value of Q(s). 
According to (241, the series B(t) given by (26 1 converges in the circle |f| < 1/a, because B {N) ~ (-a) N N b "~ b ° as 
N —» oo. The nearest singularity of the series is located on the negative real half-axis at the point t = -1/a. Then 
the integration contour over t e [0, +oo) intersects the boundary of the circle of convergence for expression (25) at 
the point 1/a. The problem of analytical continuation of ( [26) beyond the convergence domain \t\ < 1/a can be solved 
either by the method of the conformal mapping of the complex plane or by Pade approximation method H24H . Below 
we will use the conformal mapping method, because it is controlled by HOA. Furthermore, it leads to more accurate 
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results than other methods (see [25]). In our case the position of the B(f)-function poles depends on the position of 
the invariant coupling constants in the (gi, go) plane. 

For example, let us consider a system with r — 4. There are two kinds of instantons in the model. For 
instanton containing one non-zero block we get a(l) = 3(2g t + g of) 14. Otherwise, instanton has two non-zero blocks 
and a( 2) = 3(4g\ + gf)/8- Therefore in stability sector <[9]) there are two regions in the plane (g \, go) where series for 
B(t) have different analytical properties: 

Region I: if the invariant coupling constants satisfy the condition 8gi + 3go > 0, then \a( 1 )| > \a(2)\. The nearest 
singularity of the B(t ) is t - -l/a(l); 

Region II: if the invariant coupling constants satisfy the condition 8gi + 3g2 < 0, then |a(2)| > \a( 1 j|. In this case 
the nearest singularity of the B(t ) is located on the positive real half-axis at the point t — l/|a(2)|. 

Thus, the plane (g i, g 2 ) is divided by the line 8gi + 3g2 = 0. Above this boundary the analytical properties of 
B(t)-fu actions are determined by one non-zero block instanton, under the boundary only the two non-zero blocks 
instanton influences the properties of function B(t). 

The initial values of the invariant coupling constants are located in the region I. Let us apply conformal mapping 
method for the invariant coupling constants located in this region. Usually the conformal map of the complex plane is 
chosen in the form (23l [24l l 


u(s) = 


Vl + as - 1 
Vl + as + 1 


o 


s(m) = 


4 m 

a(u - l) 2 


(27) 


The series ([25]) can be rewritten in terms of the variable u as 


B(e) = Z £(A0 g /V _ z U (N) u N , 

N>0 N >0 

n ( 28 ) 

U (0) = B {0 \ U m = B (m) (4/ a) m C^~m -1 > N > 1, 

m= 1 


then the conformal Borel map of the quantity Q looks as follows 


6(e) = 



dt t b °e 'ii{st) N . 


(29) 


Usually, the parameter bo is chosen to weaken the singularity of the Borel transform (26 1 at the point 
fixed by the relation bo = b n +3/2 li23l [24j. 

In the region II the singularity of the function B(t) is located on the positive real half-axis, thus 
mapping method can not be used. 


/ = -1/a. It is 
the conformal 


3.2. Numerical analysis of the RG-equations. 


Combining the RG equations (15 ) and resummation formula (291, we have resummed RG equations for the in¬ 
variant coupling constants 


K 

d{gi = ~gi + Y U< i N) f dt tb ° e ~' u ( £t ) N . 

w=o X 


(30) 


£;l f =o _ S'- 


Note that U {N) and u(t) are functions of the variables g,. This system of equations (30) can be solved by standard 
finite-difference method. 

The results of numerical solutions of the system ( f30| for r - 4 are shown as an example in Fig[2] at s = 1 and 
Fig[3]at e = 2. Fig. [2]shows that in three dimensional model the invariant charges trajectories starting with different 
initial values cross the boundary of the action stability domain at some value fo of the parameter Similar behavior 
is observed for different values r > 4. 
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Figure 2: Trajectories of the running coupling constants at D = 3 and r = 4; dashed line - the boundary of applicability of the resummation 
method. 



Figure 3: Trajectories of the running coupling constants at D = 2 and r = 4; dashed line - the boundary of applicability of the resummation method. 



Figure 4: The solutions of the RG equations(D = 3, r = 4) at different numbers of calculated loops: 1-loop - dotted line, 5-loops - solid line. 



Fig|4]shows how the trajectories of invariant charges depend on the order of loops calculations for D - 3. We can 
state that five-loops approximation is sufficient to ensure the loss of the action stability and accurate calculation of 4 : o. 
Moreover, numerical analysis shows that solutions of the Cauchy problem ( [30] ) are stable under small perturbations 
of initial conditions. 

It is interesting to note, that we have found the IR stable fixed point of RG equation. According to mm there 
is no IR stable fixed point for /3-functions in one, two and three-loops approximations. The fixed point appear in four- 
loops. But the five-loops corrections essentially change the position of this fixed point, so we can guarantee neither 
existence nor position of this fixed point. 

In D — 2 (e = 2) case the IR stable fixed point is found in four- and five- loop approximation too. But in difference 
to three-loops ones m only rare trajectories of the invariant charges cross the line of the action stability (|9]) according 
to Fig.Q. These trajectories are connected with very small initial values of renormalized coupling constants. 

Our calculations are valid if the invariant coupling constants are in the region I. In the region II the series can not 
be resummed by the Borel method. Finally note that a( 1 j gives the largest contribution to the HOA in the stability 
domain and in a neighborhood of the stability boundary for any r > 2. One can assume that the phase transition occurs 
near the boundary of stability. For this reason resummation process can be made only for a(l). 


4. The phase transition description. 

The loss of the action stability is usually considered as a mark of the first-order phase transition. But obviously, 
it is not possible to claim that £o defines the first-order phase transition temperature; only metastable states appear in 
the system at £ = To answer the question when the new state of the system (with a condensate) in fact becomes 
stable, i.e., to determine the phase transition temperature, more accurate analysis is necessary. 

Because the interaction terms (~ x 4 ) of the action (10 1 are not positively defined now, we have to take into account 
the next term (~ x 6 ) of the “bubble” expansion of the action { 5 ]). 

Let us consider an effective action (10 1 with an additional F 2 = ti(x^x) 3 term. In renormalization procedure in 
4 - s scheme F? will be considered as a composite operator of canonical dimension A 3 = 6 - 3e. Also, there are 
composite operators F 2 = tr(xX) 2 U '(XX) and F 1 = (tr ^ 1 ^) 3 with the same canonical dimension as F 3 , therefore they 
may be mixed in the process of renormalization. Thus the term AojFj/36 must be included in the effective action, here 
Toy are bare homogeneous sources. One can define the set of renormalized parameters A, using Aoj = Zj^A^M 2 ^ 2 for 
such extended model, Z ; /, = 6 jk + Zj ; /s + ()(1 /e 2 ). Matrix Z is a function of the variables g;. Similar to the (13 1 we 
can write RG functions for A< 


Paj = -(2s 


Q 

• 2)Aj + Aigk- — Z]j. 

ogk 


( 31 ) 


The matrix Z was calculated in the framework of three-loop approximation, see appendix. One-loop approxima¬ 
tion of matrix Z leads to the following results 


At, = 2(1 - S)A\ +g 1 


^li(r 2 — r -H 14) + A 2 (r - 1) + 


+ ^82 [4i(r - 1) + A 2 ], 


(32) 


Pa 2 - 2(1 - s)A 2 + -g\ [/l 3 ( 6 r - 9) + A 2 (r~ - r + 38)] + -g 2 [6 A\ + A 2 (r - 2) + 3 /I 3 ], (33) 

(34) 


Pa 3 = 2(1 - e)/l 3 + ~^~giA 2 + -g 2 [A 2 (r - 4) + 4A 2 ], 


the rescaling of charges gj —> gi /1 6tt 2 is assumed. 

Let us mark that the full family of the composite operators with the same canonical dimension in the D — 4 
dimensional space must be taken into account for an accurate calculation of the F, operators renormalization. 

This family also includes the operators 

h = tr(A* + A*), / 41 = tr(Ax + xx + x) + tr(x + AxX + X), M = tr(d,x + xd,X + X) + tr(; Y + diXX + diX), ( 35 ) 

/ 43 = tr(Ax + x) tr(x + X) + trCr + A^) trCr + ^), f 44 = tr (d,x + x) tr(d,X + x) + tr( X + d>X ) tt(x + d<x). Us = K(d i x + d i x) K(x + X) 
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Figure 5: Thermodynamics potential as a function of order parameter: a - disorder state, b - metastable state, c — “superfluid” state. 


in addition to the F, operatorsd. The canonical dimensions of these operators are d[f 2 \ - D + 2, c/[/ 4 ,] = 2D - 2. But 
in the analysis presented we limit ourselves by the considerations of the F, operators only. The contributions of the 
operators ([35} will be discussed below. 

It was shown in liT4l that p = (x) is an order parameter of phase transition in the model considered. A non-zero 
value of (3 leads to phase transition to the superfluid phase. 

The value for magnitude /3 can be calculated by minimization of free energy —T. In the framework of the Landau 
mean field theory this functional can be written in the form 


- r = T tr/j’/J + 8 -f (tr/T ) 2 + 8 -f tr ^ (upp f f + ^ tr (pp f f tr/?/? f + ^ tr (/3j8 f f . (36) 

Schematically it can be represented in the figure (Jsj). The variables /3 and //' have the Pfaffian’s form (To} . For the 
extrema conditions at the phase transition point we get 

^-r = o, -^-—r = o, r = o, v; = o,...,r/2. (37) 

d/3j dp* J 

Obviously, the loop corrections to equation ([36} contains IR singularities. This singularities can be taken into account 
using RG method. This procedure leads to the fact that charges goj, Aoj in ( [36} must now be replaced by the invariant 
charges gj, Aj, which in turn depend on the parameter r. After such processing, the contributions of higher loops give 
only e-corrections to the mean field theory results. Let us introduce Zj = fij/ M d '‘ , s = t/M 2 , where dp = 1 — e/2 is 
canonical dimension of the field p. Then RG equations for the invariant variables are 


a - 

% * i= 2 + r / 

&J 9 

(38) 

^-2 + r / 

^i = 0 — 

(39) 

d * Zj ~ ~ Zj 2 + y r ’ 

• Zj\ e= o ~ z i ’ 

(40) 


Finally, if we combine previous equations with condition ([37} and (36 1 , we get 


\u 2 = ~ 


2 «gi + gi 


2 4n 2 Ai + 2nAz + /I 3 ’ 


(41) 


9 (2»g! + g 2 ) 2 

16 4n 2 7li + 2nA2 + T 3 ’ 
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T = 


(42) 










Figure 6 : Trajectories of “effective” coupling constants A = 4n 2 Ai + 2nA 2 + /I 3 (dashed line) and G = 2ngi + g 2 (solid line) at D = 3 and n = 2. 



Figure 7: Trajectories of “effective” coupling constants A = 4irA\ + 2nA 2 + /1 3 (dashed line) and G = 2ng\ + g 2 (solid line) at 


U = 1 and n ■■ 


n is number of non-zero blocks. Thus, as r decreases, the invariant charges intersect the boundary of the stability 
domain and new solution ( 141 ) of stationary equations ( | 37 ) appears. This phase have two non-zero blocks, n = 2 . 
Equation (42 1 determines the transition temperature r,. In order to solve equation ( | 42 ) it is necessary to know solutions 
of RG equations (38 1 and (39 1 . As before, the RG equations must be resummed. Similar to © we can rewrite the 
equations (39 1 


deb = 


'M^o - Ai ' 


2 e - 2- 


A i + J]s N AjLfCg u -g 2 ), 


N =0 


The HOA for iJ.' 1 are needed for the Borel resummation too. Our analysis in Sect.2 shows that calculation of L' 

coefficients is connected with the renormalization of six-point 1PI Green functions (~ ( ViV) 6 ) which include one 
insertion of composite operators F,- ~ ( ViV) 6 , hence L (Ar> /( ViV) 6+6 ~ B (N> /( ^fN) 4 . Indices structure is irrelevant for 
the HOA. " 

Thus, we can resum the RG equations ( |43) by the formula (29 1 . The results of numerical computations are shown 
in figures[6]and[7] This allows us to solve the equation (42 1 . As the result the root of this equation differs only a little 
from fo obtained in Sect.3. Remember that fo demonstrates a weak dependence from initial values of the coupling 
constants g 


(43) 


(AO 


Let us discus here a possible contributions of / operators (35 1 to the results obtained. There are some reasons why 
we have not calculated these counterterms. 

First, we can state that these contributions are relatively small compared with F because / operators are more IR 
irrelevant in the real space dimensions D — 2,3 then F according to the canonical dimensions mentioned above. Then 
the corresponding invariant charges will be oppressed by the first terms in the RG equations similar to ( |43) for these 
variables. 

Second, it can be simply shown by the instanton analysis presented that the high-order contributions of / operators 
are small in 1 /N compare with these of F. 

And third, the calculations of the renormalization of the full family of composite operators F and / is rather 
technically difficult now. To calculate the full renormalization constants matrix up to e 3 corrections one needs to 
consider six-loop diagrams. It is not worth while to start these calculations, as our results show that the first order 
phase transition takes place and the influence of F operators on it’s temperature is rather small. 
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Then we can state that in the model considered the first-order phase transition takes place at a temperature higher 
than the predictions for continuous phase transitions. 

To estimate the temperature difference in D - 3 case we have calculated numerically that £o = T t /g\ ~ 2 4 3 in a 
wide range of values gi ~ 10 s 4- 0.1, here r, is the root of the equation (42 1 . Natural to assume that the chargesare of 
the same order of magnitude gi ~ g 2 « 1. Then the renormalization constants Z,i have the form Z,i - \ + O(gj). In 
this approximation the ratio Z r /Z^ equals to 1. This leads to the relation 


T o/ g ()2 = ft- 


(44) 


The integrals over momenta and sum over frequencies a> s (|7][8]) can be reduced to the one-dimensional integrals. But 
the RG-approach used in our article give us an opportunities to calculate different values for small r only. Then it is 
sufficient to calculate these parameters 00 using the approximation /36 » 1. It can be found in this approximation 


g 02 


7 vfP 

8 (nT) 2 


f(3). 



co 


7 vfP 2 f P 

96(nT m ) 2 


f(3), 


(45) 


with corrections ~ 0(([35) '). Here vp = mpp I(2 tt 2 ) is 3Z>-density of states at the Fermi level, pp is the Fermi 
momentum. 

Near the transition point tq can be estimated as 


f3vp AT 
2cq Tq 


(46) 


where To is the continuous phase transition temperature determent by the usual approach JT|. Combining (44 45 46 1 
we get estimation for the temperature difference between the first order phase transition and Tq 


AT _ y 69 1 27T 6 / T \ 4 
7f(3) \7>) ' 


(47) 


5. Conclusions 

In contrast to the case of the electron systems (r — 2, r - number of spin degrees of freedom) where continuous 
phase transition takes place, our investigation has shown that in systems with high spin fermions (r > 4) critical 
fluctuations destroy stability of the system (see Fig.2). In such systems the first order phase transitions take place in 
space dimension D = 3. These results were obtained by means of renormalization group analysis with e-expansion 
up to the fifth-loop order of perturbation theory and subsequent Borel resummation. It should be noted that five loop 
calculations are indispensable to be sure that the first order phase transition takes place. 

The temperature of transition to the superconducting or superfluid phase was estimated for the systems under 
consideration. Three loop RG analysis for composite operators, which are similar to 1 ) 3 in the Landau-Ginzburg 
functional, was performed for estimation of this temperature. It was revealed that transition temperature is higher 
than the theoretical estimation based on the continuous phase transition formalism for the same model. The obtained 
difference in temperatures is rather small (see expr. (|47j)). But it should be kept in mind that approach used in the 
present work is applicable for the small deviations from the phase transition temperature only. Thus, in either case, 
we can guarantee that the difference in the phase transition temperature is not lower than our estimation. 

As for 2D systems one can state that the five loop approximation is not sufficient to determine neither the phase 
transition type nor the phase transition temperature. The last is an excellent example in favor of further development 
of the high-loop calculations. 
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6. Appendix 

The results obtained for coupling constants renormalization in three-loop approximations are presented here. Five- 
loop results for /3 gj , y T and three-loop results for/T are available by kalagovg@gmail.com. 

1 3 

Pg, = ~eg l + 4 O 2 ~r + 8 )g\ + (r- 1 )gig 2 + -g\- 

- 2(3r 2 - 3r + 14)g 3 - ^(r - 1 )g 2 g 2 - 2(5r 2 - 15r + 92 ) gl £ - |(r - 2 )&+ 

+ ^ [33r 4 - 66 r 3 + (955 + 480<T(3))r 2 - (480^(3) + 922)r + 2960 + 2112^(3)] g\+ 

+ [79r 3 - 158r 2 + (1397 + 768<T(3))r - 1318 - 768^(3)] l? 3 g 2 + (48) 

+ ^ [3 r 4 - 12r 3 + (576^(3) + 3355)r 2 - (1728^(3) + 7568)r + 9216^(3) + 14788] g]gl+ 

+ ^ [60r 3 - 32 lr 2 + (2943 + 1152£(3))r - 2304^(3) - 4092] t?1 g? + 

+ [(96^(3) + 193)r 2 - (576^(3) + 891)r + 1536^(3) + I860] g\, 

1 , 

Pgi = ~eg2 + ^( 2 r - 5)gj + 3 £ig 2 - 

- ^(r 2 -lr + 2 (% 3 - i(ll r - 20 )g 2 2 g l - 2(5 r 2 - 5r + 82)g 2 g 2 + 

+ — l — [26r 3 - (383 + 96f(3))r 2 + (2459 + 1152£(3))r - 4060 - 2688^(3)] g\+ 

1024 ( 49 ) 

+ ^ [(96^(3) + 182)r 2 - (963 + 576<T(3))r + 1536^(3) + 1937]g 3 g 1+ 

+ yyt [ _70r3 + 1 lr2 + ( 6423 + 4608<T(3))r - 8064(73) - 10366] g]gl+ 

+ [- 13r4 + 26r 3 + (192^(3) + 355)r 2 - (368 + 192f(3))r + 3284 + 2688^(3)] g 3 g 2 

Tr = ~j(r 2 ~r + 2)gi - l(r - l)g 2 + 2 ( r 2 - r + 2)g 2 + -(r - l)gig 2 + 2 ( r 2 - 3r + 4)g 2 

- -2(15r 4 - 3Or 3 + 156r 2 - 141r + 222)g 3 - -2(15r 3 - 30r 2 + 126r - 111 )^ 2 (50) 

256 128 

- -2-(r 4 - 4r 3 + 403r 2 - 960r + 888 )gig? - 2—(90r 3 ~ 321r2 + 987r “ 888 )^ 

J. . 1 J. \J Z* 1 
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